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Abstract 

We revisit some topics of classical thermostatistics from the perspective of 
the nonextensive optimal Lagrange multipliers (OLM), a recently introduced 
technique for dealing with the maximization of Tsallis' information measure. 
It is shown that Equipartition and Virial theorems can be reproduced by Tsal- 
lis' nonextensive formalism independently of the value of the nonextensivity 
index. 
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I. INTRODUCTION 



Tsallis' thermostatistics (TT) [|T]-|5|] is by now recognized as a new paradigm for statistical 
mechanical considerations. It revolves around the concept of Tsallis' information measure S q , 
a generalization of Shannon's one, that depends upon a real index q and becomes Shannon's 
measure for the particular value q = 1. 

The Equipartition theorem, first formulated by Boltzmann in 1871, and the Virial one, 
due to Clausius (1870) || are two pillars of classical physics. They were discussed within 
the TT framework in J7|, in terms of the Curado-Tsallis nonextensive normalizing treatment 
0, and also in |§, by recourse to escort probabilities 0. In the present work we are going 
to revisit these classical subjects in connection with a recently advanced scheme for dealing 
with Tsallis' thermostatistics flQ]]Tl|| , that seems to yield illuminating insights into classical 
themes. 

Tsallis' thermostatistics involves extremization of Tsallis' entropy by recourse to the 
celebrated technique of Lagrange. A key TT ingredient is the particular way in which 
expectation values are computed || and, in such a respect, several proposals have been 
considered during the last ten years ||. No matter which recipe one chooses (for q ^ 1), 
classical phenomena are always reproduced by TT in the limit g -> 1 0. Recently, a new 
algorithm (to be henceforth referred to as the OLM-one) has been advanced to such an 



extremizing end that diagonalizes the Hessian associated to the Lagrange procedure flO 
and yields the optimal Lagrange multipliers associated to the input expectation values. A 
diagonal Hessian enormously facilitates ascertaining just what kind of extreme the Lagrange 
method is leading to |10| . 



A key point in TT considerations is the following one: the entropy constant k is usually 
identified with Boltzmann's ks- However, the only certified fact one can be sure of is 
"fc — > ks for q — > 1" fl3|, which entails that there is room to choose k = k(q) in any suitable 



way. It is seen |L0| that if one chooses 

k = k B Zf\ (1) 



where Z q stands for the partition function, in conjunction with the OLM formalism ||10|| . 
the classical harmonic oscillator determines a specific heat C q = ks, so that the classical 
Gibbs' result arises without the need to invoke the limit q — > 1. Additionally, this treatment 
is able i) to reproduce thermodynamics' zero-th law [[LI]], a feat that had eluded previous 
TT practitioners, and ii) the classical g-independent value for the mean energy of the ideal 



gas [pT2| . It is then not unreasonable to conjecture that many other classical results might 
be obtained by TT without invoking the q — > 1 limit. 

In the present effort we re-discuss the equipartition theorem and the virial one from an 
OLM viewpoint. We also revisit the ideal gas problem as an application. More precisely, we 
perform an OLM treatment of the subject following the canonical ensemble strictures. 

II. BRIEF REVIEW OF OLM FORMALISM IN A CLASSICAL SCENARIO 

The most general classical treatment requires consideration of the probability density 
p(x) that maximizes Tsallis' entropy [|I|,p|,|l3"f 

S q _ 1 - [dx p 9 (x) 

fc~ q-1 ' {) 

by recourse to the Lagrange's technique, subject to the foreknowledge of M generalized 
expectation values 

(in \\ - / rfx P g (x)Oj(x) , v 

WW*- /dxp*(x) ' [6) 

where Oj(x) (j = 1, . . . ,M) denote the M relevant dynamical quantities (the observation 
level [|14j]), q € 3? is Tsallis' nonextensivity index, x is a phase space element (N particles in 
a .D-dimensional space), and k is the entropy constant, akin to the famous Boltzmann one 
ks, employed in the orthodox statistics. 

Tsallis' normalized probability distribution 0, is obtained by following the well known 
MaxEnt route ]T5| . Instead of effecting the variational treatment of involving M + 1 



Lagrange multipliers Xj (associated to constraints given by the normalization condition 



3 



together with the M equations @), the OLM technique follows an alternative path flO 



with Lagrange multipliers A^: one maximizes Tsallis' generalized entropy S q (|2]) fi],|2|.|T3| 
subject to the modified constraints ("centered" generalized expectation values) [|I|,|i~0H 

J dx p(x) -1=0 (4) 
JdxpWfaW-UOj)),) =0. (5) 



The resulting probability distribution reads |10 



PW = — 7? , (6) 



where 



Z g =/rfx/(x) — . (7) 



and 



M 



/(x) = 1 - (1 - g) £ A} (0,-(x) - (( 0j )) q ) (8) 

3 

is the so-called configurational characteristic. 

Although the Tsallis-Mendes-Plastino (TMP) procedure originally devised in over- 
comes most of the problems posed by the old, unnormalized way of evaluating Tsallis' gen- 
eralized mean values fOlJTtl, it yields probability distributions that are self-referential, which 
entails some numerical difficulties. The complementary OLM treatment of jT(j] surmounts 



these hardships. Inspection of (§) shows that the self-reference problem has vanished. 
One shows in |10] that the relation 



J rfx p*(x) = Z]~\ (9) 

valid under TMP, still holds. Eq. (||) allows one to connect the Lagrange multipliers Xj of 
the TMP procedure with the corresponding OLM Aj via [Kj 



and to write the entropy as JT^ 



S q = k \n q Z q} 



(11) 



where ln g x = (1 — x l q )/(q — 1) has been used. 
If now, following IJ we define 



In, Z' q = \n q Z q - A; <(0^ 



(12) 



and additionally use |TU 



k' = k Z\-\ 



(13) 



then 10 



dS n 



d((O j )) ( 



k'K 







r(KZ' q ) =-((Oj)) q . 



(14) 
(15) 



Equations ( |T4|) and (|l^) constitute the basic Information Theory relations to build up 
Statistical Mechanics a la Jaynes j0|. Notice that, when k is given by (|l|), Equation ([13]) 
leads to k! = ks 0. 

Remembering that || 



d((O j )) l 



kX 



j- 



(16) 



one is straightforwardly led to [10 



k'X'j = kXj, 



(17) 



which entails that the intensive variables are identical in both alternative pictures, TMP 
and OLM ]TJ . 

As a special instance of Eq. (|TTD, for the Canonical Ensemble it adopts the appearance 



k>{3> = k{3 = -. 



(18) 



Looking at Eq. (|T8| ) one gathers that the temperature T is the same for both sets of 
Lagrange multipliers. 

Before tackling the equipartition theorem, some preliminary results are needed, that we 
discuss next. 



III. NORMALIZATION CONSIDERATIONS 

Replacing (|6]) into @ we obtain for the partition function the expression 

Z q = fdxf{x)&, (19) 
and, comparing it to (^), we deduce the (at first sight surprising) relationship 

D = Jdx /(x)t^ - J dx /(x)^« = 0, (20) 
valid for all q. That this is indeed so can also be gathered by first recasting ( ^0|) in the form 

V = |dx[/(x)^(l-/(x))], (21) 
and then using Eq. (H) to obtain 

V = (1 - 9) E*5 / rfx /^ - <(^»J > ( 22 ) 



that, remembering (||), has to vanish identically on account of (|J). 
We introduce now the factor 

m - ™i = i, (23) 

/dx/(x)i-« 
in order to be in a position to use it later on. 

It is to be stressed that F(q) = 1 only within the so-called "normalized" framework of 
||. If one uses, instead, the Curado-Tsallis normalization ||, this is not so (F(q) ^ 1). The 
OLM formalism does employ the normalized treatment. It must be pointed out that the 
F(q) = 1 -result has also been obtained within a purely quantal, Green function scheme by 
Lenzi, Mendes, and Rajagopal fl7j . 
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IV. GENERALIZED EQUIPARTITION AND VIRIAL THEOREMS 



In classical statistical physics, a Hamiltonian dynamical system is described by an ap- 
propriate phase space probability distribution p(ri,Pi). Just one assumption will be made 
on the probability density (PD): that it depend on the phase space variables {ri,pi) only 
through the Hamiltonian H(ri,pi). 

Tsallis' normalized probability distribution || is obtained by maximizing Tsallis' gener- 
alized entropy S q given by Eq. (0), subject to the constraints flUJlOf 

dnp( ri , Pl ) = i (24) 

J dn P \ ri , Pl ) (H{r hPi ) - {(H)) q ) = 0, (25) 

where the generalized expectation values 

//m\ ! dVLpijr^p^Hjr^p,) 

m) * = jdnp^) (26) 

are assumed to be a priori known. They constitute the macroscopic information at hand 
concerning the system. dQ stands for the corresponding phase space volume element 



DN 



dfl = (l/(N\h DN )) J] dridpi, (i = 1 . . . DN), 



(27) 



i=i 



where h is the linear dimension (i.e. the size) of the elementary cell in phase space, and we 
assume / T\i=l d r i = V N with V the system's volume. 



The resulting probability distribution reads JT0| 



pin, Pi ) = Z- 1 [l - (1 - q)(3' (H(r hPi ) - ((H)) ( 



1-9 



(28) 



where 



dVL 



l-(l- g )/3' (H(n, Pi )-((H)) ( 



1-9 



(29) 



Let A(q it pi) denote a generic dynamical quantity. Its generalized mean value is given by 

jdSiAinjpt) \i-0d-q) (h-((h)) c 



((A)) c 



9 
1-9 



JdH l-p(l-q)(H-((H)) 



9 

1-9 



(30) 



Let us now assume that our Hamiltonian H can be split into two pieces according to 



H = g + g 



(31) 



where g is a homogeneous function of degree 7 of L canonical variables of which, say u, are 
generalized coordinates while the remaining ones are p = L — v generalized momenta 



9 = 9(r 1 ,...,r u ,p 1 ,...,p li ), 



(32) 



and go does not depend upon these variables. According to Euler's theorem we have 



19 = J2 r i( d 9/ dr i) + J2Pi( d 9/ d Pi) 

i=l i=l 



(33) 



so that the generalized mean value of g reads 

1 



09)) q = £ EiHdg/dn))), + HW*,))), 

/ U=l i=l 

We shall now discuss in some detail one generic term of this equation. Consider 

„ ra /a uv Jdnr h (dg/dr k )[l-l?(l-q)(H-((H)) t 
((r k (dg/dr k ))) g = - 



(34) 
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(35) 



fdn [l-p(l-q) (H-{{H)) q) 

a mult i- dimensional (2DN) integral. Let us evaluate the integral over r k ranging between 
r a and r&. These values are given by the well-known Tsallis' cut-off condition [Q: the 
probability distribution vanishes in those regions of phase space that would make ( p8|) a 
negative quantity. For the numerator of ( P5|) we have 



J 



dr k r k (dg/dr k ) 1 -/?'(!- g) (F - ((H)) 



1 

l-q 



(36) 



so that 



((r k {dg/dr k ))) c 



J ...J J dr 1 ... dr k _idr k+ i ...dp 



DN 
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(37) 



fdn [l-(3'(l-q) (H-((H)) q 

In order to obtain J we proceed to an integration by parts. To do so we first notice that, 
on account of Equations ([H]) and (|32"D 



S 



_d_ 

dr k 



l-/3'(l-g) (H-((H)) t 



l-q 



(38) 



which, since the integrated part will vanish because of the above mentioned cut-off condition, 
leads to 



l-0{l-q) (H -((&)) 



l-q 



(39) 



Insertion of fl3*9|) into fl3~T|) yields 

1 Jdtt \l-/3'(l-q) (H-((H)) a 
((r k (dg/dr k ))) c [ 1 ^ ^ " 
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JdQ l-L3'{l-q)[H-({H)) 



i i 

l-q 



(40) 



but, on account of (|23D one has 



((r k {dg/dr k ))) q 



1 



It is apparent that each term in the sums appearing in 
the type fl4"l|) , so that 



((g)) q = A = -k B T 



(41) 

will yield a contribution of 

(42) 



which is a generalized version of the equipartition theorem ||. We have thus arrived to 
a new an interesting result. Contrary to current belief, we see that the classical result is 
attained, independently of the q-value. 

According to the canonical equations of motion, dH/dri = —pi. Hence ( |4"T| ) leads to the 
statement 



t \i=l 



-DNk B T 



(43) 



which is Clausius' virial theorem . Again, no dependence upon q is to be detected. 

As a simple application, consider now the classical ideal gas. We necessarily reproduce 
the classical results, by virtue of the above considerations: we deal with N particles confined 
within a D— dimensional box of volume V. Thermodynamical equilibrium at temperature 
T is assumed. The Hamiltonian is 
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N 2 N 
j=l i=l 

where U wa ii is a contribution due to the constraints (walls of the container). The constraint 
forces come into existence when the gas particles collide with the walls. It thus follows that 

E*V^ = -E^-Fi:L (44) 

i=l UL i i=l 

where F^ all stands for the force on the i-particle due to the walls of the box. 

Following now a well-known argument M, it is possible then to express the last member 
of the above equation in terms of the volume V and pressure P 

N 

-E*V F il, = -W r-(-ds) = P vds (45) 

j_j J wall J wall 

= P y (V • r)d 3 r = DPV, 
where ds denotes a (vector) surface element of the box in Eq. (H5). We get 



N QTT 

^■i~ = DPV (46) 

Now, by taking the canonical average and using the virial theorem (unmodified by the 
nonextensive scenario) fl43|), we obtain 

PV = Nk B T, (47) 

i.e., the equation of state for the perfect gas. No g-dependence is detected. 

For the ideal gas @,|,0U], the total energy E is a homogeneous quadratic function 



(7 = 2) of DN momenta, which allows one to write, according to (|42|) 



((H)) q = \DNk B T. (48) 



V. CONCLUSIONS 



Classical thermostatistics has been the subject of the present effort. 
We have tackled some key issues, namely, 
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• Virial theorem, 

• Equipartition theorem, 

• Equation of state of the ideal gas, 

and shown that they are, contrary to present belief, reproduced by Tsallis' thermostatistics 
independently of the value adopted by the index q. The present work lends further credence 
to the hypothesis that most important classical statistical results might be reproduced by 
Tsallis' thermostatistics for all g-values. 
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